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ABSTRACT 

We report the discovery of a partial ~ 2° diameter non-thermal radio shell coincident with 
Taz, the pulsar wind nebula (PWN) in the error box of the apparently variable 7-ray source 
3EG J1809— 2328. We propose that this radio shell is a newly identified supernova remnant 
(SNR G7.5— 1.7) associated with the PWN. The SNR surrounds an amorphous region of ther- 
mal X-rays detected in archival ROSAT and ASCA observations putting this system in the 
mixed-morphology class of supernova remnants. G7.5— 1.7 is the fifth such supernova remnant 
coincident with a bright GeV source, and the fourth containing a pulsar wind nebulae. 

Subject headings: gamma rays: observations — supernova remnants — pulsars: individual(Taz) — radio 
continuum: ISM — X-rays: ISM — ISM: individual(G7.5-1.7) 



1. Introduction 

The possibility of supernova remnants (SNRs) 
being sources of observable GeV emission has long 
been suggested by t heoretical models o f cosmic- 
ray acceleration (see Baring et al. 1999, and ref- 
erences therein) in SNR shocks and by statis- 
tical associations with ca t aloged 7— ray sources 



( Sturner fe Dermerl [l995 lEsposito et all [l996; 



Yadigaroglu fe Romanil 19971 ). However, the di 
rect linking of a particular SNR with a specific 
7-ray source detected by the EGRET instrument 
on board the Co mpton Gamma-Ray Observatory 
has b een elusive (|Dohertv et al . 2003; iTorres et al 



2003) (note we only refer to the results of the su- 
pernova blast as SNR, not isolated pulsar wind 
nebulae (PWN) which are sometimes referred to 
as "plerionic" or "Crab-like" SNR). The EGRET 
sources with the smallest and most robust error 
boxes are not positionally coincident with regions 
of SNRs which show direct evidence of high-energy 



particle acceleration through synchrotron emis- 
sion of X-rays. Although there are five previously 
known radio SNRs that contain within their shells 
unidentified EGRET sources which are bright 
abov e 1 GeV (as defined by lLamb fe Macombl 
19970 . in no case can hard X-ray emission from 
the SNR blast wave be convincingly asso ciated 
with the 7-ray source ([Roberts et al.ll2001al here- 
after, RRK). 

Interestingly, four of these SNRs ar e of the sub- 
class of mixed- morphology SNR fSNR lRho fe Petre 
19981 ) which have centrally concentrated thermal 
X-ray emission within a radio shell. Three of 
these (CTA 1, IC 443, and W44) also contain 
young, energetic pulsars which are actively pro- 
ducing high-en ergy particles as evidenced by their 
wind nebulae dBrazier et al. 19981 Halpern et al 



20041: lOlbert et al.ll200lUHarrus et al.lll996f) . It is 
therefore unclear as to whether the GeV emission 
is produced at the SNR forward shock, the interior 
where the thermal X-ray emission comes from, or 
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has to do with the pulsar or its wind. 

While the source of GeV emission remains un- 
clear, at TeV energies a growing number of sources 
coincident with PWN and youn g shell SNR have 
b een discovered by H.E.S.S. (jAharonian et al 



2006cl) . Several of the TeV emitting PWN are 
associ ated with EGRET sour ces, notably the 
Crab dWeekes et alJll989h . Vela (lAharonian et al 



2006al ). and t he two nebulae in the Ko okaburra ra- 
dio complex (|Aharonian et al.ll2006bl ). Two clear 
detections of s hell SNRs w ith H.E.S.S. have also 
been reported ( Funk 2006). In these cases, the 
TeV emission has been shown to be associated 
with regions of non-thermal X-rays in the shell. 
A few H.E.S.S. sources coinciden t with composite 
SNR have been repo r ted a s well dAharonian et al 



2005; iBrogan et ail 120051 : iHelfand et all 120071) . 



but in these cases the TeV emission appears to 
be associated with non-thermal X-rays from the 
central PWN, not the shell. Two of the mixed- 
morphology SNR associated with bright GeV 
sources, IC443 and W28, have recently been de- 
tected as sources of TeV emission, but it is unclear 
what, if any, connection there is between the GeV 
source and the TeV source ( Albert et al.l 2007; 



H. E. S. S. Collaboration: G. Rowell et alJl2007h . 

Here we report the detection of a non-thermal 
radio shell with the Very Large Array ( VLA) at 
90 cm surrounding the PWN associated with the 
EGRET source G eV J1809— 2327 /3EG J1809 -2328. 
Arch ival ROSAT (jTrumperlll982[) and ASCA dTanaka et 
1994) data show a region of diffuse thermal X-ray 
emission between the radio shell and the PWN 
suggesting this source is a new mixed-morphology 
SNR. 

2. Analysis of Archival X-ray Data in the 
Region surrounding the Taz PWN 



3EG J1809-2328 is one of the b rightest sources 
of > 1 GeV emission in the Ga laxy ( Hartman et al 
19991: lLamb fc Macomblll997i) . It has a well de- 
termined position (Fig [lj and apparently vari- 
able emi ssion above 100 M eV on timescales of 
months (jNolan et al.1 l2003h . ASCA imaging 
of the EGRET error box revealed an extended 



hard X-ray source (RRK), which (jBraie et al 



l2002h subsequently resolved into an apparent 
rapidly-moving PWN (RPWN) and a young stel- 
lar cluster using a short Chandra observation 



fsee lKaspi et~ai1l2006t iGaensler fc Slanejl2006L for 
general reviews of PWN). Radio imaging with 
the VLA and observations with VMM-Newton 
confirm the RPWN nature of the X-r ay source 
((Roberts et al.ll2001bl : iBraie et al.ll2002L Roberts 
et al. in preparation). Due to the radio nebula's 
distinctive funnel shape (presumably imposed by 
its motion), its unusually powerful 7-ray emission, 
and the growing tradition of naming PWN after 
animals, this nebula is sometimes referred to as 
"Taz" (short for tasmanian devil). 

The X-ray and radio morphology of Taz, if in- 
terpreted as being created by a bowshock, suggests 
a birthsite to the northwest. A ROSAT PSPC 
exposure corrected 0.1-2.4 keV image obtained 
through Sky ViewfH shows a large region of soft X- 
ray emission in this direction which we will refer 
to as G7.5-1.7 (Fig. [J). This inspired a 38 ks 
observation with the ASCA satellite on 1996-09- 
24 that has not previously been published. Fol- 
lowing the imaging and spectral fitting procedures 
outlined in RRK, we produced combined particle 
background subtracted, exposure corrected images 
ofthc3EG J1809-2328 and G7.5-1.7 fields from 
the ASCA GIS data in the 0.7-2 keV, 2-10 keV, 
and 4-10 keV bands (Fig. [J). The higher energy 
bands show a fairly bright, compact source of hard 
X-rays in the region of soft thermal X-ray emission 
which is not evident in the ROSAT image. The 
peak of this source is at RA = 18 ,i 08 m 56.3 s , Dec 
g -23°09'57" (J2000) in the GIS 2-10 keV image; 
ASCA images have a nomi nal positional error of 
~ 24" (jGotthelf et al.ll2000h . We will refer to this 
hard source as AX J1808.9-2310. 

We extracted spectra from two regions in the 
field of the GIS instruments, one containing the 
hard source and one containing the soft dif- 
fuse emission. Since there is emission essentially 
throughout the ASCA field, we used a background 
region from the neighboring field which contains 
the Taz PWN. We fit the soft extended emission 
to a thermal model, choosing an absorbed vnei 
non-equilibrium model from XSPEC to facilitate 
comparison with other Mixed- Morphology rem- 
nants, at first keeping the abundances fixed to 
solar. To help constrain the soft end of the ther- 
mal emission spectrum, we simultaneously fit the 



1 The NASA Sk yView Virtual Observatory is located at 
http: / /skyview.gsfc. nasa.gov/ 
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3EG J 1809-2328 




Right Ascension (J2000) 

Fig. 1. — Multiwavelength view of G7.5-1.7. MSX 8.3 
/an (20" resolution) is in red, ROSAT PSPC 0.1-2.4 keV 
X-rays (smoothed with an 70" FWHM Gaussian) are in 
blue, VLA 90 cm radio (168" X 120" resolution) is in green. 
The contours are the 68%, 95% and 99% uncertainty on 
the position of 3EG J1809-2328. The circles show X-ray 
spectral extraction regions (see fig. 2). The dotted cyan 
lines delineate the nominal shell of G7.5— 1.7 while the 
dashed red line shows the approximate outline of the Lynds 
227 dark cloud, whose edge is bright in the infrared image. 
The 90 cm VLA image has been corrected for primary beam 
attenuation. 



ROSAT PSPC spectrum extracted from the the 
same region. We obtained an adequate fit us- 
ing a single temperature with solar abundances. 
The result is a temperature kT = 0.61+"'^ keV 
(90% confidence region), an ionization timescale 
t = 1.54+04? x 10 11 scm -3 and an absorption 
n H = 1.5±o'.l x 1021 cm" 2 . Note that the 
ROSAT and ASCA spectra seemed to differ 
somewhat near the upper end of the ROSAT pass- 
band (~ 2 keV) and so these values are some- 
what dependent upon which energy bins from 
ROSAT were included in the fits. For exam- 
ple, ignoring the ROSAT data above 1.6 keV 
resulted in an improved chi-squared and a some- 
what lower absorption of ~ 9 x 10 20 cm~ 2 . Fitting 
the ASCA data by themselves resulted in consis- 
tent values, and allowing the abundances to vary 
did not result in a significantly better fit. We also 
tried the ray and vray equilibrium thermal mod- 
els, but these resulted in statistically inferior fits. 
The spectral values obtained are typical of temper- 
atures and timescales fit from ASCA observations 



of ot her mixed- morphology SNR ( Kawasaki et al 



2005). Since the thermal X-ray emission covers 



an area much larger than the field of view of the 
ASCA GIS and it is somewhat unclear what is the 
total extent even in the ROSAT images, an ac- 
curate total thermal flux measurement is difficult. 
However, using the fluxed ROSAT maps and ex- 
trapolating from our fit spectra, we estimate the 
total unabsorbed bolometric thermal flux from the 
remnant to be on the order of 10~ 10 ergcm~ 2 s^ 1 . 

We then fit the hard source (AX J1808. 9-2310) 
to an absorbed power law plus an absorbed vnei 
model, fixing the latter to the values above al- 
lowing only the normalization to vary. The re- 
sult is an absorption nn = 1.9+04 x 10 22 cm -2 , a 
spectral index V = 2.57±oj? and a 2-10 keV nux 
F x = 3.9 ±0.1 x 10~ 12 ergs cm" 2 s -1 . Since there 
is no evidence for a source of additional foreground 
absorption in the mid-infrared or optical (see be- 
low) , the order of magnitude higher absorption to- 
ward this source suggests it is an unrelated back- 
ground object and/or contains significant intrinsic 
absorption. The total Galactic nu in this direc- 
tion as estimated us ing the HE AS ARC nn t ool 
and the H I maps of Dickey fc Lockman ( 1990l ) is 
- 6 x 10 21 cur 2 . 
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3. Radio Observations 

We observed the region around Taz at 90 cm 
with the VLA on May 28, 2004 in the DnC config- 
uration (project code AR 547) for 4.5 hours (on- 
source). We used 32 channels over a bandwidth 
of 12.6 MHz in LL and RR polarizations. Us- 
ing AIPS and standard wide-field l ow frequency 



imag ing techniques (see for example iBrogan et al 
20041 and references therein), the resulting im- 



age has an effective central observing frequency of 
324.84 MHz, an effective bandwidth of 9.38 MHz, 
a beam size of 168" x 120" (P.A.=72°) and an RMS 
level of about 20 mJy/beam. The VLA FWHP 
primary beam at 325 GHz is 2.5°. Due to spatial 
filtering caused by the missing short spacings, this 
image is not sensitive to smooth structures larger 
than - 1°. 

The 90 cm image shows a partial shell of 
emission surrounding the thermal X-ray emis- 
sion (Fig. [TJ. Comparison with the 8.3// MSX 
image of the region shows no infrared excess as- 
sociated with this shell, sug gesting a non-thermal 
origin of the radio e mission iCohen fc Greenll200ll 



Brogan et al. I l2006h . We convolved the 90 cm im- 
age with a 2-D Gaussian to produce a 258" res- 
olution image for comparison to the Bonn 11cm 



Survey image for this region ( Reich et al. 1984h . 



We made two spectral index maps from the 90 cm 
and 11 cm images (assuming S"), one using the 
standard 11 cm background subtracted image, 
and one using an 11 cm image that was passed 
through a "high-pass" filter in order to mimic the 
large-scale spatial filtering of the 90 cm VLA im- 
age. In both spectral index maps, the Taz PWN 
and the new shell both have significantly steeper 
spectra than known thermal sources in the region. 
Overall, the spectral index map created from the 
filtered 11cm Bonn data (Fig. [3j) is steeper, as ex- 
pected, and provides a better match to the spectra 
of other known SNR in the field of view. The spec- 
tral index of the shell ranges from a ~ —0.3 to 
—0.7, with an estimated uncertainty of ±0.3 at 
any particular position. The large uncertainty is 
dominated by the inherent uncertainty in mea- 
suring flux from a large, diffuse object with an 
interferometer, and not the noise in the images. 

The center of the radio shell is at roughly 
I = 7.54°, b = -1.90° near the northern edge of 
the Taz radio nebula, and the radius is r ~ 0.82°. 




Right Ascension (J2000) 

Fig. 2. — ASCA GIS X-ray image of G7.5-1.7. Red is 
0.7-2 keV, green is 2-10 keV, blue is 4-10 keV (~ 1' resolu- 
tion) . The circles show the extraction regions for the X-ray 
spectral fits of the supernova remnant, the hard source, and 
the background region. 




t" r 

-22°oo' ^._o sjs 3& / SNR 

20' 
40' 
-23°00' 
20' 
40' 
-24°00' 
20' 

1 B n 1 4 m 12 m 10 m 08 m 06 m 04 m 
J2000 Right Ascension 

Fig. 3. — Spectral index (assuming S") map of 
G7.5— 1.7 and the surrounding region derived from an 
11 cm Bonn image and the 90 cm VLA image. Before cal- 
culating the spectral index, the Bonn image was passed 
through a "high-pass" filter to mimic the spatial cover- 
age of the VLA interferometric data (which are not sen- 
sitive to smooth strucures larger than 1°). The 90 cm 
image was first convolved to the 258" resolution of the 
Bonn image. The input 11 and 90cm images were masked 
at 30 and 120 mjy beam -1 , respectively. Contours from 
the 168" X 120" resolution 90 cm image are superposed, 
the contour levels are 60x(l, 2, 3, 4, 5) mjy beam -1 . 
Several previously id entified H II regions and SNRs are 
labeled for referenc e (Griffith ct al. 1994; Lockma n et alj 
ll996l : lBrogan et al.ll2006h . 
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The shell has a somewhat smaller radius of curva- 
ture (r ~ 0.45° center at I ~ 7.51, b ~ -1.42) on 
the side nearer the Galactic plane, centered in the 
middle of the thermal X-ray emission. Assuming 
the origin of the expanding shell lies somewhere 
in between these two centers and near the symme- 
try axis of Taz, we will refer to the radio shell as 
G7.5-1.7. 

We have also obtained high resolution (3.6" x 
3.2") 20 cm VLA data of a much smaller field of 
view centered on the Taz nebula which will be 
discussed in detail in a future paper that concen- 
trates on the PWN (Roberts et al. in prep.). For 
the purpose of the present paper it is interest- 
ing to note that there is a 13 mJy 20 cm radio 
source at RA(J2000) = 18 h 08 m 57.59 s , Dec(J2000) 
= -23°09'45.4" that is coincident with the hard 
ASCA source AX J1808.9-2310. The radio 
source is elongated and well fit by a 7.4" x 1.4" 
Gaussian, suggestive of a barely resolved double 
jet source. 

4. Discussion 



G7.5— 1.7 has the classic properti es of a mixed- 
morp hology supernova remnant ( Rho &: Petrel 
1998). The thermal X-ray emission is concen- 
trated between Taz and the Galactic plane side 
of the radio shell. In both the Bonn 11 cm and 
our VLA 90 cm images, between 1/2 and 3/4 of 
a radio shell surrounds Taz. The side towards the 
Galactic plane is brightest, while the side away 
from the plane is too faint to be sure of any struc- 
tures, suggesting a strong density gradient in this 
direction. Comparing to the MSX 8.3/im image, 
the radio shell is located within a region of low 
diffuse mid-infrared emission. Numerical models 
of the evolution of mixed-morphology remnants 
suggest they are a common stage in the growth 
of remnants at ages of ^10,000-100,000 yr if they 
are expanding in h igh-density regions of the ISM 
( Tillev et al. 2006). In this case where there is 
an apparent density gradient, the thermal X-ray 
emission is concentrated near the bright side of 
the shell, towards the Galactic plane, rather than 
around the apparent center of the shell. 

The position of the center of G7.5— 1.7 suggests 
it is the birth remnant of Taz. While the overall X- 
ray absorption of Taz is significantly higher (nn ~ 
1.8 x 10 22 cm~ 2 , RRK) it is also moving into or 



behind the Lynds 227 dark cloud. XMM-Newton 
imaging shows some parts of the X-ray nebula ar e 



absorbed more than others ( Roberts et al. 2006h . 



and a preliminary spectral analysis suggests a 
range of nn ~ 0.7 — 2.0 x 10 22 cm~ 2 . If we as- 
sume Taz and G7.5— 1.7 are associated, then the 



distance to Lynds 227 of 1.7 kpc (|Oka et al.lll999h 



can be considered a lower limit to G7.5— 1.7. 

From the above estimates of the center of the 
radio shell, and the pulsar wind nebula line of 
symmetry, we estimate the birthsite of Taz to 
be I ~ 7.53°, b 1.68°. The average trans- 
verse space velocity is then implied to be vt ~ 
200 (di, 7/^50) km s , fairly typical for an isolated 
pulsar. Here di.7 is the distance in units of 1.7 kpc 
and £50 is the age in units of 50 kyr. The distance 
from this estimated birthsite to the pulsar is about 
half the average distance to the shell (Rj, ~ 0.62°), 
which is roughly where a bowshock is expected to 
begin forming aroun d a pulsar moving within a de - 
celerating SNR shell (jvan der Swaluw et al. I l2004h . 

If we accept the association with the pulsar and 
hence the distance to Lynds 227 as a lower limit 
on the SNR distance, and we further assume it has 
been undergoing Sedov expansion for the majority 
of its existence (i.e. R oc t 0A ) then we can estimate 
the current expansion velocity of the shell to be 
v s ~ 150(rfi. 7 /t 50 ) km s _1 . Note that if the SNR 
were considerably younger (eg. t < 15,000 yr), 
or considerably further away (d > 5 kpc), then 
the velocity of the radio shell would be similar 
to very young SNRs. For example, only a few 
degrees away on the sky is the X-ray and radio 
bright shell of the ~ 2000yr old Type II remnant 
SNR Gil. 2— 0.3 which has a well measured ex- 
pansion vel ocity of ~ 700 km s" 1 at a distance 
of ~ 5kpc ([Roberts et all 120031: iTam fc Robert^ 
2003). For strong shocks, the temperature be- 
hind the shock, which in the case of Gil. 2— 0.3 
is kT ~ 0.6 keV, is proportional to the square of 
the shock velocity. Although the flux from such a 
shock depends somewhat on the density and com- 
position of the ISM as well as the equilibrium state 
of the electrons with the ions, such a hot shock 
would be difficult to hide from ROSAT given the 
rather bright radio emission and overall low ab- 
sorption. However, at the estimated velocity of 
~ 150km/s given our nominal distance and age, 
the bulk of the shock emission would be emitted 
well below 1 keV where absorption is more impor- 
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tant and ROSAT is much less sensitive. Therefore, 
the lack of bright X-ray emission associated with 
the radio shell is consistent with the interpretation 
of G7.5— 1.7 being at a distance of ~ 2 kpc and 
an age similar to other mixed-morphology SNR. 

G7.5— 1.7 is now the fifth known mixed- 
morphology supernova remnant that is coinci- 
dent with the error box of one of the 25 bright- 
est sources of GeV emission in our Galaxy, and 
the fourth of these to contain a pulsar wind neb- 
ula. Two of these SNRs also contain known 
TeV sources, and several other PWN associ- 
ated with bright G eV sources have TeV emission 



(jFunk et al.l 120071 ) , so it is quite plausible that 
Taz/G7.5 — 1.7 also has TeV emission. However, 
in no currently known case is the GeV emission 
convincingly arising from the same spatial loca- 
tion as the TeV emission. This suggests that GeV 
and TeV emission may have separate sites of ac- 
celeration or represent particle populations with 
very different ages. 

As is the case with the PWN in CTA1 and 
W44, Taz is actually contained within the GeV 
error box, while the radio shells are outside the 
error box. In all three of these cases, there is 
some evidence of moderate variability of the 7-ray 
emission, with the sourc e containing Taz s howing 



the strongest evidence ([Nolan et al 1 l2003h . The 



W44/PSR B1853+01 system is very similar in that 
there is clearly an RPWN seen in both radio and 
X-rays at the edge of the thermal X-ray emitting 
regio n. It also has similar X-ray ((Kawasaki et al 



2005) and 7-ray properties (RRK). The rapid mo- 
tion of the pulsars through the clumpy dense in- 
terior of these remnants provides a natural vari- 
ability mechanism for emission coming from the 
immediate environment of the pulsar. However, 
PSR B1853+01 is only moderately energetic with 
a spin-down power E — 4.3 x 10 35 erg s" 1 . There- 
fore, if this is where the variable 7-rays from W44 
are produced, then it suggests the interaction of a 
pulsar wind with the unique hot and dense envi- 
ronment of the interior of mixed-morphology rem- 
nants may allow for very efficient 7-ray produc- 
tion. 
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